Beyond the rainbow: effects from pion back-coupling 
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We investigate hadronic unquenching effects in light quarks and mesons. To this end we take 
into account the back-couphng of the pion onto the quark propagator within the non-perturbative 
continuum framework of Schwinger-Dyson equations (SDE) and Bethe-Salpeter equations (BSE). 
We improve on a previous approach by explicitly solving both the coupled system of DSEs and BSEs 
in the complex plane and the normalisation problem for Bethe-Salpeter kernels depending on the 
total momentum of the meson. As a result of our study we find considerable unquenching effects in 
the spectrum of light pseudoscalar, vector and axial-vector mesons. 
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I. INTRODUCTION 

In hard scattering events mesons and baryons can 
be viewed as bound states built up from partonic con- 
stituents, i.e. quarks and gluons. This picture changes 
at low energies, where hadronic effects play a more promi- 
nent role in the nonperturbative structure of hadrons. Of 
particular importance are pion cloud effects which e.g. 
have a direct impact on the spin structure of the proton 
Thus they need to be incorporated in bound-state 
calculations aiming at a realistic description of mesons 
and baryons. 

For the spectrum of light pseudoscalar mesons, it is the 
axial-vector Ward-Takahashi identity that governs the 
pattern of dynamical symmetry breaking, rather than 
the intricacies of confinement or infrared specifics of 
QCD. Within the framework of Schwingcr-Dyson (SDE) 
and Bethe-Salpeter (BSE) equations, this realisation 
has led to much phenomenological success in describing 
pseudoscalar and vector mesons [5, 's', '2], with bound- 
state masses, decay constants and electromagnetic form- 
factors determined 0, Q with surprising agreement to 
experiment f?\. 

The success of this description partly comes from the 
vector-vector coupling of the rainbow-ladder truncation, 
where the resultant spin-coupling to these 5- wave mesons 
is dominant. For P-wave states however, such as the 
scalar and axial-vector mesons, a richer tensor structure 
is required from the quark-gluon vertex to properly ac- 
count for the observed spectrum. To this end, the scalar 
part of the quark-gluon vertex, whose importance regard- 
ing confinement has been highlighted in [8|, may be the 
most relevant for introducing a spin-dependence to the 
interaction. We expect any beyond-the-rainbow study 
that generates such terms to have an effect on the spec- 
trum of light mesons. 

However, there arise further complications in describ- 
ing the scalar and axial- vector sector. Not only are mod- 
els employing the rainbow-ladder truncation too attrac- 
tive, yielding masses of 800-900 MeV 0, Q, but are these 



mesons even qq bound states? For the scalars we have 
tetraquarks, multi-meson molecules, glueballs and the in- 
evitable hybrids as candidates (see e.g. Indeed, 
these exotics are not restricted to the scalar sector, for 
the ai ma y b e well described as a coupled-channel meson 
molecule 113, amongst other equally plausible descrip- 
tions [13, Il4| . Clearly there is much work to be done 
in providing a dynamical partonic description of these 
'mesons', not least in understanding the fundamental in- 
teractions between quarks and gluons. 

As a further step in this direction we concentrate in this 
work on improving our model description of qq bound- 
states. Based on an approximation scheme for the quark- 
gluon interaction developed in [l3| and we study the 
pion back-reaction on the quarks and light mesons in 
the covariant Schwinger-Dyson/Bethe-Salpeter approach 
to Landau gauge QCD [13, [l^- As a consequence we 
take into account pion cloud effects in the light me- 
son spectrum. To this end we solve the coupled sys- 
tem of quark Schwinger-Dyson equations and Bethe- 
Salpeter equations for pseudoscalar, vector and axial- 
vector mesons in the complex plane. This provides for 
an important technical progress compared to the calcu- 
lations in Ref. [lB| , where the so-called real- value approx- 
imation has been used. As a further technical complica- 
tion we then have to solve the normalisation problem for 
Bethe-Salpeter amplitudes generated by a kernel depend- 
ing on the total momentum of the meson. 

The paper is organized as follows. In Sec.|TTl we sum- 
marise our approximation scheme for the quark-gluon in- 
teraction developed in detail in [H, [l^ . We specify our 
model interaction and discuss the technical complications 
that appear in the normalisation condition for the bound- 
state amplitudes. We explicitly verify that our interac- 
tion satisfies the axial-vector Ward-Takahashi identity. 
In Sec. mil we then provide the details of our calculation 
and present our results. We conclude with an outlook in 
Sec. IIVI Some technical discussions are relegated to an 
appendix. 
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II. THE QUARK-GLUON INTERACTION 
A. General truncation 

In the past, a widely used practical approximation 
scheme for the coupled system of quark-SDE and meson- 
BSE is the rainbow-ladder approximation. Its most im- 
portant property is that it satisfies the non-singlet axial- 
vector Ward-Takahashi identity (axWTI) which guaran- 
tees that the pion in the chiral limit is both a Gold- 
stone boson and a bound state of massive constituents. 
This scheme has a history of remarkable successes (sum- 
marised e.g. in However, there are also shortcom- 
ings that limit the credibility of such an approximation. 
Consequently, several efforts have been made to extend 
this scheme, see e.g. [H, IM El, [13, HI Hi ■ In the fol- 
lowing, we will continue the strategy employed in [l^ to 
extend the rainbow-ladder scheme by taking into account 
pion cloud effects in the quark-gluon interaction. These 
effects can be accounted for by a particular class of in- 
teraction diagrams that only appear in the unqucnchcd 
theory. Our strategy is complementary to the one fol- 
lowed in (25} . where a different class of unquenching ef- 
fects have been investigated. The aim there was to pro- 
vide for hadronic intermediate states in bound-state cal- 
culations which generate a finite width of meson spectral 
functions. It is also complementary to the investigations 
reported in [l^, [2^, where unquenching effects in the 
gluon polarization have been considered. 

We start with the SDE for the quark propagator, given 
diagrammatically in Fig. [TJ With the fully dressed in- 




FIG. 1: The Schwinger-Dyson equation for the fuUy dressed 
quark propagator 

verse quark propagator S~^{p) — i-^ A(p^) + B{p'^) and 
its bare counterpart Sq^{p) = ij/ + m, the equation is 
given by 

S-\p) - Z2S^\p) 

with k = p - q, the Casimir Cp = {N^ -l)/{2Nc) 
and the renormalization factors Zip = Z2/Z3 of the 
quark-gluon vertex, Z2 of the quark propagator and Z3 
of the ghost dressing function. The quark dressing func- 
tions A{p'^) and B{p^) can be recombined into the quark 
mass M{p'^) — B{p'^)/A{p'^) and the quark wave function 
Zf{p'^) = 1/A{p'^). The quark self-energy depends on the 
fully dressed quark-gluon vertex T^{q,k) and the gluon 




FIG. 2: The approximated Schwinger-Dyson equation for the 
quark propagator with effective one-gluon exchange and one- 
pion exchange. 



FIG. 3: The Bethe-Salpeter equation corresponding to the 
quark self-energy of Fig. [2] The up-down arrow indicates an 
averaging procedure of the pion-exchange diagram with its 
counterpart where the upper pion-quark vertex is bare and 
the lower one dressed. 



propagator 

, , / kakA Z(k^) 

D,Ak)^(S,.~^j^, (2) 

with the gluon dressing function Z{k'^). Throughout this 
paper we will work in the Landau gauge, which guaran- 
tees the transversality of ([2]). The explicit expression for 
Z(fc^) used in this work will be given in Sec. Ill B[ 

The key quantity in the SDE for the quark prop- 
agator is the quark-gluon vertex. It satisfies its own 
Schwinger-Dyson equation which contains various basi- 
cally unknown four-point functions. The structure of 
the Yang-Mills part of this equation and the resulting 
pattern of dynamical chiral symmetry breaking has been 
discussed recently in detail in Ref. Q . Here we are inter- 
ested in the hadronic contributions to this vertex. These 
have been discussed in Ref. fTs*! and [16*1, where two 
slightly different approximations in terms of one-pion ex- 
change have been constructed. Whereas the scheme of 
[Tsl resulted in a very strong back-reaction of the pion 
onto the quark propagator in disagreement with cor- 
responding lattice results, the modified scheme in [l^ 
seems to represent these effects quantitatively correctly. 
We therefore resort to the latter scheme here, whose di- 
agrammatic expressions for the quark-SDE and the cor- 
responding Bethe-Salpeter equation for quark-antiquark 
bound states are given in Figs. [5] and [3] 

The effective self-interaction of the quark is thus given 
by a pure Yang-Mills part and a hadronic part approxi- 
mated by one-pion exchange. The dressed coupling of the 
pion to the quark line is hereby described by its Bethe- 
Salpeter vertex function, whereas the bare one is given 
by T^'^^Z2 with the generators of flavour SU(2). 

Certainly it is important that the explicit form of both 
the Yang-Mills part and the pionic part of the interaction 
is such that the axial-vector Ward-Takahashi identity is 
satisfied. Naturally, as a first guess for the Yang-Mills 



3 



part one takes a rainbow-ladder form of the interaction, 
i.e. one chooses a model for the fully dressed Yang-Mills 
part of the quark-gluon vertex of the form 



(3) 



where we denote the quark momenta by pi and p2 and 
the gluon momentum by ps . The explicit expression for 
pYM^^I jg detailed below. A similar form has been used 
in [23, [231 , where quenched lattice results for the quark 
propagator have been reproduced. Note, however, that 
Q involves only the j^-part of the tensor structure of the 
full Yang- Mills part of the vertex. It has been shown in 
the analysis of the full quark-gluon vertex of Ref. [8] that 
such a model cannot capture all essentials of dynamical 
chiral symmetry breaking. Indeed, from the results of 
this work we will draw the very same conclusion from 
our calculation of the light meson spectrum. 

In general one can decompose the quark-gluon vertex 
T'^{pi,P2,P3) into 12 different tensor structures, given by 



r,.(pi,P2,P3) 



12 

E 

i=l 



K{Pl,P2,P3)Ll{Pl,P2,P3) ■ (4) 



The details of this basis are given in Ref. [29|. Here, 
we only note that — 7j^; the explicit forms of the 
other structures are not needed. In principle, all twelve 
tensor structures are important in the infrared and in- 
termediate momentum regime |30| . We have verified 
by projection methods that the pion-exchange diagram 
in the vertex-SDE indeed contributes to all these struc- 
tures. Therefore the resulting one-pion contribution in 
the propagator SDE Fig. [2] generates unquenching effects 
in all components of the quark-gluon interaction. 



An explicit expression for the one pion-exchange ker- 
nel that satisfies the axial-vector Ward-Takahashi iden- 
tity has been constructed in [15]. Here we give a similar 
construction which takes into account the modifications 
suggested in [l6| . The resulting SDE for the quark prop- 
agator is given by 



S-\p)^Z2So\p) 



Z(fc^)FYM(fc^ 

A:2 



(27r)4 



Z2l5S{q)TT, ( ^^-^■,q 



Z2l5S{q)T^ [ ^^-^■,p- q 



(5) 



with k = p — q and D-,r{k) = (fc^ + M^)^^ being the pion 
propagator. The factor 3 in front of the pion contribution 
stems from the flavour trace and represents contributions 
from 7r+ , 7r_ and ttq . These are treated on equal footing 



in the isospin-symmetric limit adopted here. 

The corresponding expression for the Bethe-Salpeter 
equation for light mesons reads 



r^iP^ P)=j [Kj^rsiP. k; P) + Kt^AP. k; P)] [s{k+)T^^\k-, P)S{k^) 



with the kernels K^J^g and Kf^°^g given by 



K^i,{q.p;P) 

Kf!::^s(<i,p;p) 



g^Z{P)T^^{k^)Z^F 



[ri] 
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Ts{p;P) 
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Ta{p;P) 
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TABLE L The four Dirac structures describing a meson of 
spin J = 0. For the pseudoscalar, these have an associated 
factor of 75 to account for parity. Circles indicate that the 
component is part of the basis decomposition of the meson, 
with fiUed circles indicating that the component is multiplied 
by {p- P). This ensures that for equal-mass constituents a 
Chebyshev decomposition of even order need only been em- 
ployed for the scalar amplitudes Fi. 



general basis for each desired meson amplitude, 



-,1+ 
(9) 

where the components T^'^'' are given in Tables [J and |lTl 
In particular, the pion is given by the following form 



(p; P) =75 [Fi [p- P)-irF2 (p; P) 

{p ■ P) Fsip; P) - ,y ] F^ip; P) 



(10) 
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TABLE II: The eight Dirac structures describing a meson of 
spin J = 1. For the axial- vectors, there is an additional factor 
of 75 to account for parity. Circles indicate that the compo- 
nent is included in the basis decomposition of the meson, with 
filled circles indicating the component is multiplied by (p ■ P). 
This ensures that an even Chebyshev decomposition is suffi- 
cient for equal-mass constituents. The subscript T indicates 
transversality with respect to the total momentum. 



B. Model details 

Next we need to specify the details of our model in- 
teraction. For the Yang-Mills part of the interaction 
we choose two different model ansatze that have been 
employed already in previous works. Here, the gluon 
dressing function Z{k'^) from and the Yang-Mills part 
r^^(fc^) of the quark-gluon vertex Eq. ^ are combined 
and represented by a single function. For the Maris- 
Tandy (MT) model |^] this function is given by 

_^ 27r7m 

log (^+(i + 97A|cd)') 

X [l-cxp(-gV[4™2])]'j ,(11) 



Here T'^^^\p] P) is the Bethe-Salpeter vertex function of a 
quark-antiquark bound state specified below, whereas F:^ 
with flavour index j is the corresponding vertex function 
of a pion. The momenta fc+ — k + P/2 and k- = k — P/2 
are such that the total momentum P of the meson is 
given by P = fc+ — fc_ and the relative momentum k — 
(fc+ -I- fc-)/2. The Latin indices {t,u,r,s) of the kernels 
refer to colour, flavour and Dirac structure. 

The Bethe-Salpeter vertex function F(p;P)(''^ can be 
decomposed into at most eight Lorentz and Dirac struc- 
tures. The structure is constrained by the transforma- 
tion properties under CPT of the meson we wish to de- 
scribe 31|. In particular, our pseudoscalar, scalar and 
vector have quantum numbers of 0~, 0^ and 1~, re- 
spectively. The axial-vector can be parameterised in two 
different ways, dependent on its transformation under 
charge conjugation, J^*-^ = 1++ and 1^ . 

Taking this into consideration, we can write down a 



with 



mt = 0.5 GeV 
7,„ = 12/(33- 



r = 
Aqcd 



e"-l 

= 0.234 GeV. 



This model has been employed successfully in various cal- 
culations of meson spectra, form factors and decay con- 
stant within a rainbow-ladder truncation; see (TtI . [2^ for 
overviews of the results. The remaining parameters D, 
Lj and the quark mass are fitted to pion observables, and 
detailed in Sec. IIIII 

The second model we employ for the Yang-Mills part of 
the quark-gluon interaction represents the infrared struc- 
ture of both the gluon dressing function Z{k'^) and the 
quark-gluon vertex F"*'^(fc^), in quenched approximation 
(up to effects due to hairpin diagrams) as calculated in 
Refs. 0, Isij . It has been used to show that infrared diver- 
gences in vertex functions can account for the topological 
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mass of the Ty'-meson [33|, and is given by 



z = 



I / afit (fc ) \ 



+ ^^QCD 



r™(A;2) = 



-1/2-K 



(12) 



fc2 



d2+fc2 d2^(fc2_rf2)' AqCD+fc' 

-25 . 

47r / 1 



A2 

^^QCD 



^ V"qcd/ 



P _ A2 

^^QCD 



with the gluon momentum k^, the one-loop value 7 ~ 
{-ISNc + 4Nf)/{22Nc - -iNf) for the anomalous dimen- 
sion of the re-summed gluon propagator, the correspond- 



ing one 6 



for the ghost dressing function and 



We 



- 0.2 at the renormalisation scale /i^ = 170GeV^. 
use Aq(^]3 — 0.52 GeV^ similar to the scale ob- 
tained in Ref. ^32i] . The infrared exponent k has been 
determined analytically in [s^, and is given by k = 
(93 - VI2OI) /98 ~ 0.595. Note that this form is slightly 
modified to that reported in [16. . .33] , allowing more free- 
dom in separating the scales in the quark-gluon vertex 
and the gluon propagator. The parameter di sets the 
strength of the interaction in the infrared. The parame- 
ters d2 and da determine the location and size of a term 
added to give additional integrated strength in the in- 
termediate momentum regions. These three parameters, 
together with the current quark mass, are constrained 
by pion observables and the generation of a topological 
charge in the pseudoscalar isosinglet channel. The run- 
ning coupling from the ghost-gluon vertex, aatip^), is pa- 
rameterised such that the numerical results for Euclidean 
scales are accurately reproduced (32| : 



An 



-p2/A2 



QCD 



QCD 
1 



(13) 



ln(p2/A|cD) 



p2/A2 



QCD 



1 



Here /?o — (H^c ~ 2A'^y:)/3, and as{0) is the fixed point 
in the infrared, calculated to be as{0) — 8.915/A^c for 
our choice of k. In this work, we will take Nf = 2. 

The resulting two effective interactions from the Maris- 
Tandy model (fTTj) and our soft-divergence model (fT2|) are 
plotted against each other in Fig. U) The ultraviolet run- 
ning of both interactions are the same except for differ- 
ences in the scale employed. Whereas the Maris- Tandy 
model used a scale corresponding to the MS scheme, we 
are working in a momentum subtraction scheme where 
the scale is larger by more than a factor of 3. Both in- 
teractions share an enhancement in the intermediate mo- 
mentum region which provides the necessary interaction 
strength for dynamical chiral symmetry breaking. The 



I Maris-Tandy 



10 10' 10 10 10 10" 10 

[GeV"l 

FIG. 4: A log-log plot of the effective quark-gluon interaction 
Zr^'^' for the two models we employ. Differences in the UV 
behaviour are attributed to the different scale Aqcd in use. 
Note both exhibit enhancement in the intermediate momen- 
tum region. 



change of slope in the infrared region of our interaction 
is related to the top ological mass component of the 77' in 
the chiral limit [33l |. 

As for the pion part of our interaction, we approxi- 
mate the full pion Bethe-Salpeter wave function PH]) in 
the quark-SDE (O and the kernel of the BSE p by the 
leading amplitude in the chiral limit given by [2| 



ri(p;P)=T^75 



(14) 



Here B(jP) is the scalar dressing function of the quark 
propagator in the chiral limit. This approximation omits 
the back-coupling effects of the three sub-leading am- 
plitudes. From the calculation of ref. [3], where full 
back-coupling has been evaluated in a real value ap- 
proximation, we expect this omission to lead to an er- 
ror of only a few percent as concerns meson masses and 
of about 10-20% as concerns decay constants. Consid- 
ering the exploratory character of our work this seems 
certainly tolerable. The replacement of the leading phys- 
ical pion amplitude Fi{p;P) at — 138 MeV by the 

(normalised) chiral limit value Fi {p; P) sa ^^f~^ ^^^^ 
correct on the few percent level as demonstrated in the 
appendix. Note that we use this approximation only for 
the internal pion which mediates the interaction. The 
great advantage of the approximation (|14l) compared to 
the full back-coupling performed in Ref. J^] is that we 
can then fully take into account the quark propagator 
in the complex plane as necessary input into the Bethe- 
Salpeter equation. 



C. Decay constants and normalisation 

The Bethe-Salpeter amplitudes for the pseudoscalar 
mesons can be used to obtain the corresponding weak 
decay constants /^r, which will prove useful in fitting pa- 
rameters to observables. In order to calculate these one 
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must firstly normalise the amplitudes. The normalisation 
condition is derived by demanding that the residue of the 
pole in the four-point quark-antiquark Green's function 
(from which the BSE is derived) be unity [s^ls^]- Using 
conventions such that — 93 MeV, and for momenta 
partitioned equally between the two constituents, it reads 

^^^■-2^tr/|^ (15) 

3 (^r;(fc, -Q)S{k + F/2)ri(fc, Q)S{k - P/2) j 

where = — is fixed to the on-shcU meson mass, 
the trace is over Dirac matrices and the Bethe-Salpeter 
wave-function x is defined by 

Xi(fc; P) = S{k + F/2)ri(fc, P)S{k - P/2) . (16) 

The conjugate vertex function T is given by 

T{p,~p)^cr^{~p,-p)C'\ (17) 

with the charge conjugation matrix C = —7274- The 
leptonic decay constant characterising the pion coupling 
to the point axial field is subsequently given by [bt! ] 

3 f k 

h = Z2j^tvj ^r.(fc,-P)5(fc+)75f 5(fc_) , 

(18) 

where again the trace is over Dirac matrices, and k^ = 
k+P/2, k- = k — P/2. There exist analogous expressions 
for the vector mesons where an additional factor of 
1/3 must be included due to summation over the polari- 
sation tensor. 

One may write a similar equation to (jlSp . which cor- 
responds to the residue of the pseudoscalar vertex: 

r^= Z2Z„,iivj .0^r^(k,~P) S{k+)j5S{k^) . (19) 

The axWTI imposes a relationship between these two 
residues, known as the generalised Gell-Mann-Oakes- 
Renner relation which must hold at and beyond the 
chiral limit: 

Uml = (m„(^^) -|- TOd(M^)) , (20) 

where m„ , rrid are the masses of the up and down quarks 
at the renormalisation point /i^, and in this work con- 
sidered to be degenerate. Confirming this relation serves 
as a check of both our numerics and indicates how well 
our kernel satisfies the axWTI. We find excellent agree- 
ment with errors smaller than 1% for pion masses up to 
~ 600 MeV. 



III. NUMERICAL METHODS AND RESULTS 

A. Solving the BSE 

We solve the BSE for the meson amplitude F^^^ (p; P) 
via matrix methods using the following procedure. First, 
we project our BSE onto the scalar amplitudes of our 
meson decomposition. This gives rise to either four or 
eight coupled integral equations for the Fi {p; P) . To 
make manifest the angular dependence of the amplitude 
functions, we treat the total momentum P^ as a param- 
eter and expand the function as a series of Chebyshev 
polynomials in the angle p ■ P = p ■ P/\pP\ 

F^P) = Y.i^^Ftip'^P'mf^) . (21) 
fc 

The functions F,^{p'^;P'^) are projected out through use 
of the orthonormal properties of the Chebyshev poly- 
nomials. With the angular dependence made explicit, 
we can evaluate numerically the two non-trivial angles 
appearing in the integration measure. We cast the re- 
maining radial integral in the form of a matrix equation 
by matching the external momenta to the radial loop 
momenta, p^ = at the abscissae of our integration 
nodes. Thus the amplitude of our BS equation is pro- 
jected onto the decomposition F^{p'j; P^). Schematically 
we are solving 

F = AK-F, (22) 

for the column vector F as a parametric equation in 
P^ — —A'P with M the mass of the meson. A bound- 
state corresponds to solutions with A = 1, which as an 
eigenvalue problem is equivalent to satisfying the condi- 
tion det (1 — K) = 0. Because of the tractable nature of 
our unquenching prescription, the computational effort 
required to determine the Bethe-Salpeter amplitude in- 
creases from a few seconds to a few minutes. However, 
we see in the next section that substantial effort is re- 
quired to calculate the normalisation condition with any 
kernel involving momentum exchange of the meson. 

B. Calculating the normalisation 

The first term of ([15]), independent of the kernel, is 
easily evaluated since the quark propagator derivatives 
are directly calculable in our approach; either by apply- 
ing the derivative to the quark-SDE, and solving the re- 
sulting integral equation, or by finite difference methods. 
The second term is substantially more complicated for 
two reasons. First, we must determine the derivative of 
our exchange kernel with respect to the total momentum 
dependence - a tractable but algebraically complicated 
task. The second problem is the double integration over 
two four-momenta, k and q. Decomposing the measure 
into hyperspherical coordinates, where none of the angles 



7 



may be evaluated trivially, results in the need to evalu- 
ate an eight dimensional integral numerically. Resigning 
ourselves to this computationally intensive task, we eval- 
uate the integral of (fT5|) for — Q"^ ± 2"A, and employ 
the appropriate Richardson improved centered difference 
formula. Typically, we use a four-point rule for which 
n = 0, 1 though a two-point rule is sufhcient, provided 
we choose A optimally. Unfortunately each double inte- 
gral must still be evaluated, which we tackled using both 
Monte-Carlo and Gaussian quadrature techniques. For 
the latter, a small-scale computing cluster is mandatory 
for even moderate grid sizes. Monte-Carlo performs only 
marginally better in terms of necessary computing time, 
since the statistical errors must be minimised such that 
the numerical derivative is meaningful. With due care, 
however, both methods obviously yield the same result. 

In the next section, we give results for the leptonic 
decay constants for the pion and rho with and without 
consideration of the normalisation double-integral. For a 
physical pion, the difference is only about 3% and smaller 
still for the rho. However, as we reduce the pion mass and 
approach the chiral limit the necessity of including the 
full normalisation condition becomes readily apparent. 



C. Results 

The parameters of the two models we employed - that 
of Maris- Tandy and our model - were tuned such that 
pion observables were well reproduced with the inclusion 
of the hadronic exchange kernel. For our interaction, 
which is capable of generating a topological charge, the 
eta/eta' splitting provides an additional constraint. 

For the Maris- Tandy model, we find a unique choice 
of w and D for a given quark mass by fitting to TOjt 
and /tt. We attempted to maximise the rho meson mass, 
though we did not exhaust the available parameter space 
due to the complexity in evaluating the full normalisa- 
tion condition. Consequently we employ lo — 0.37 GeV 
and D = 1.45 GeV^, with a quark mass of 3.7 MeV at 

= 19 GeV. 

The parameters of our model interaction were fit in 
a similar fashion, with the additional step of calcu- 
lating the topological charge for each parameter set. 
As a result, we obtained good agreement with experi- 
ment with the choice di = 1.45 GeV^, ^2 = 0.1 GeV^ 
and ds = 3.95 GeV^. Here we choose a quark mass 
m(yu) = 2.6 MeV defined at the renormalisation point 
= 170 GeV^. 

In Table IIIIl we present a summary of meson observ- 
ables calculated for these two interactions, with and with- 
out the inclusion of the pion exchange kernel. One may 
then compare the effects of unquenching for two models 
with dissimilar behaviour in the infrared and intermedi- 
ate momentum regions (see Fig.|4]). Though the details of 
the deep infrared are generally unimportant for the pat- 
tern of dynamical symmetry breaking in light mesons, 
we expect a complicated interplay from the dynamics 





Maris- Tandy 


Our 


Model 


Experiment 




w/o pi 


incl. pi 


w/o pi 


incl. pi 


(PDG f38]) 




140 


138+ 


125 


138t 


138 




104 


93.2t 


102 


93.8+ 


92.4 






(90.2) 




(90.6) 






746 


598 


638 


485 


400 - 1200 


Mp 


821 


720 


795 


703 


776 


fp 


160 


167 


159 


162 


156 






(167) 




(165) 




Ma, 


979 


913 


941 


873 


1230 




820 


750 


879 


806 


1230 


Mr, 






493 


497 


548 








949 


963 


948 



TABLE III: BSE results for the pseudoscalar, scalar, vec- 
tor and axial-vector mesons for the Maris- Tandy and soft- 
divergence models employed with ('incl.') and without ('w/o') 
pion backreaction. The parameters marked by f are fitted to 
reproduce the pion mass and decay constant when the pion- 
exchange kernel is switched on. Consequently the model pa- 
rameters, given in the text, are modified compared to the 
original Maris- Tandy model |J] . Results for the pure rainbow- 
ladder kernel without pion back-reaction are for the same pa- 
rameter set. Decay constants given in parenthesis are those 
that would result without calculating the full normalisation 
condition of the appropriate meson. 



present in any beyond-the-rainbow truncation. 

This difference in the two interactions is seen most 
prominently for the pion. Though the effects here are 
small, as to be expected we see immediately a difference 
for the two models. Whilst for Maris- Tandy, the inclusion 
of pion-exchange is attractive giving rise to a 2 MeV shift, 
our model exhibits repulsion of the order of 13 MeV. The 
leptonic decay constants see a decrease of 10% for both 
models. The effect of including hadronic contributions is 
not then a priori known, at least for those mesons most 
sensitive to dynamical chiral symmetry breaking. 

For the remaining light mesons, the inclusion of the 
pion back-reaction generally leads to a lighter spectrum. 
In the case of the scalar meson, we see a reduction of 
150 MeV for both interactions, though the relative drop 
is significantly different at 20 and 40% for Maris- Tandy 
and our model respectively. 

The vector meson - here identified with the p - receives 
a negative shift of 100 MeV, or 12% for both interac- 
tions. In the case of Maris- Tandy, the resultant rho mass 
of ^ 720 MeV is similar to that obtained in other stud- 
ies with this interaction 4]. Here, unquenching has the 
effect of increasing the decay constant by 5%. Our in- 
teraction yields a lighter mass of around 700 MeV, with 
unquenching increasing the decay constant by just 2%. 
Since our truncation scheme does not provide a mech- 
anism for the p to decay into pions, we do not observe 
the tell-tale non-analytic behaviour of the rho mass as a 
function of the pion mass, due to the opening of a decay 
channel, as shown in Ref. [39j . 
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Next, we examine the axial- vectors oi and bi, which 
differ by their charge conjugate properties. As already 
seen, the net effect of including the pion back-reaction 
is a reduction in the meson mass. For the axial- vectors, 
we see a 70 MeV reduction in the bound-state mass for 
both charge eigenstates and in both models of the in- 
teraction. It is evident, at least with our tractable un- 
quenching truncation scheme, that unquenching effects 
do not resolve the difficulties we have in obtaining rea- 
sonable masses for these mesons. However, we note that 
the anomalous mass-splitting in calculations of the ai 
and bi in our model is half that of Maris- Tandy. In the 
course of our investigations of possible parameter sets, we 
could easily obtain near degenerate axial-vector masses 
within our interaction. 

Finally, we take a look at the masses of the 77 and 
7]'. For our model, which exhibits a soft-divergence in 
the gluon momentum, we calculate a non-zero topolog- 
ical mass following Ref. [s^. We see that the inclusion 
of the pion back-reaction on the quark propagator has 
a minimal effect on the phenomenologically determined 
mass-splittings . 

IV. CONCLUSIONS AND OUTLOOK 

Using an approximation scheme for including pion 
loop effects in the quark-gluon interaction, we investi- 
gated effects of the pion cloud in the Bethe-Salpeter 
and Schwinger-Dyson equations. By considering only the 
leading amplitude of the exchange pion, moreover taken 
from its chiral limit solution, we were able to include the 
analytic behaviour beyond the real-value approximation 
employed in the previous study . The resulting inter- 
action respects the axial- vector Ward-Takahashi identity. 

As a consequence of this technical improvement, we 
were forced to tackle the added complication of calculat- 
ing the full normalisation condition of the Bethe-Salpeter 
equation where the momentum exchange of the kernel 
need be considered. We found that this indeed gives a 
sizable contribution of the order of a few percent that 
must be properly taken into account in order to extract 
the correct leptonic decay constants. 

As a result, we found considerable pion cloud effects in 
the spectrum of light mesons. In particular we found a 
shift of about 100 MeV for the mass of the rho-meson in 
agreement with expectations from chiral extrapolations 
of quenched lattice data [ioj . 

Though this truncation scheme represents a tractable 
model for including hadronic unquenching effects, it is 



clear that spin dependent contributions from the Yang- 
Mills part of the quark-gluon vertex will have a strong 
impact upon the spectrum of light mesons. Such correc- 
tions will form a part of future investigations. 
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APPENDIX A: PION AMPLITUDE AT AND 
AWAY FROM THE CHIRAL LIMIT 

One of the key approximations here is that the leading 
Bethe-Salpeter amplitude for the pion may be well ap- 
proximated by the chiral limit result B/ f.^. It is impor- 
tant for the amplitude to feature the correct asymptotic 
behaviour at large momenta, and so for a massive pion in 
the exchange kernel we employ the scalar self-energy B 
as calculated in the chiral limit (for complex momenta). 
The differences between this and the full amplitude are 
shown in Fig. [5] and are seen to be at the level of a few 
percent. 




p [GeV] 



FIG. 5: The leading physical pion amplitude Fi{p,P) at 
rriTT ~ 138 MeV, compared to the normalised chiral limit value 
_Fi (p; P) « ^^f—^ ■ The angular dependence of this amplitude 
is negligible, and chosen here such that {p ■ P) = 0. 



[1] A. W. Thomas an d J. Lab, Prog. Part. Nucl. Phys. 61 

(2008) 219 arXivrOSOOiSf [hep-ph]]. 
[2] P. Maris, C. D. Roberts and P. C. Tandy, Phys. Lett. B 

420 (1998) 267 arXivmucl-th/9707003 . 
[3] P. Maris and C. D. Roberts, Phys. Rev. C 56 (1997) 3369 



[arXiv:nucl-th/97080"29] . 
[4] P. Maris and P. C. Tandy, Phys. Rev. C 60 (1999) 055214 

arXiv:nucl-th/9905056 . 
[5] P. Maris an d P. C. Tandy, Phys. Rev. C 62, 055204 

(2000) arXiv:nucl-th/0005015, . 



9 



M. S. Bhagwat and P. Maris, Phys. Rev. C 77 (2008) 
025203 ar Xiv:nucl-th/0612069 . [24 
J. Volmer ei al. [The Jefferson Lab F(pi) Col- 
laboration], Phys. Rev. Lett. 86, 1713 (2001) [25 
rarXiv:nucl-ex/0010009 . 

R. Alkofer, C. S. Fischer and F. J. Llanes-Estrada, Mod. [26 
Phys. Lett. A 23, 1105 (2008) arXiv:hep-ph/0607293 .; 
R. Alkofer, C. S. Fischer, ' F.^J. Llanes-Estrada [27 
and K. Schwenzer, Annals of Physics in print; 
[axXnr :0804.3042 [hep-ph]. 

R. Alkofer, P. Watson and H. Weigel, Phys. Rev. D 65 [28 

(2002) 094026 arXiv:hep-ph/0202053 . 

C. Amsler and N. A. Tornqvist, Phys. Rept. 389 (2004) 
61. [29 
M. R. Pennington, Int. J. Mod. Phys. A 21 (2006) 747 
arXiv:hep-ph/0M9265j. [30 
M. Wagner and S. Leupold, Phys. Rev. D 78 (2008) 
053001 arXiv:080L08T4 [hep-ph]]. 

J. A. Oiler and E. Oset, Nucl. Phys. A 620 
(1997) 438 [Erratum-ibid. A 652 (1999) 407] 
[arXiv:hep-ph/9702314^ . [31 
J. A. Oiler and E. Oset , Phys. Rev. D 60 (1999) 074023 [32 
"arXiv:hep-ph/9809337' . 

C. S. Fischer, D. Nickel and J. Wambach, Phys. Rev. D [33 

76 (2007) 094009 arXiv:0705.4407 [hep-ph]]^ 

C. S. Fischer, D. Nickel and R. Wilhams, ar Xiv:0807.348 6l [34 
[hep-ph]. 

P. Maris a nd C. D . Roberts, Int. J. Mod. Phys. E 12 [35 

(2003) 297 (arXiv:nucl-th/0301049 . 

C. S. Fischer,^ J. Phys. G 32, R253 (2006) [36 
arXiv:hep-ph/0605173 . [37 
C. S. Fischer, P. Watson and W. Gassing, Phys. Rev. D 
72,094025 (2005) arXiv:hep -ph/0509213 . [38 
P. Watson, W. Gassi ng and P . G. Tandy, Few Body Syst. 
35 (2004) 129 arXiv:hep-ph/0406340 . [39 
M. S. Bhagwat, A. HoU, A. Krassnigg, G. D. Roberts 
and P. G. Tandy, Phys. Rev. G 70, 035205 (2004) 
arXiv:nucl-th/0403012 . [40 
P. Maris and P. G. Tandy, Nucl. Phys. Proc. Suppl. 161, 
136 (2006) tarXiv:nucl-th/05110l7 . 
A. Bender, C. D. Roberts and L. Von Smekal, Phys. Lett. 



B 380, 7 (1996) arXiv:nucl-th/9602012\ 
H. H. Matevosyan, A. W. T homas a nd P. G. Tandy, Phys. 
Rev. G 75 (2007) 045201 arXiv:nu cl-th/0605057 . 
P. Watson and W. Gassing, Few Body Syst. 35 (2004) 99 
arXiv:hep-ph/0405287 . 

G. S. Fischer and R. Alkofer, Phys. Rev. D 67 (2003) 
094020 arXiv:hep-ph/0301094 . 

M. S. Bhagwat, M. A. Pichowsky, G. D. Roberts 
and P. G. Tandy, Phys. Rev. G 68 (2003) 015203 
[arXiv:nucl -th/0304003 . 

G. S. Fischer and M. R. Pennington, Phys. Rev. D 73 
(2006) 034029 arXiv:hep-ph/0512233 ; Eur. Phys. J. A 
31 (2007) 746 arXiv:hep-ph/0701123 . 
J. Skullerud and A. Kizilersu, JHEP 0209 (2002) 013 
arXiv:hep-pli70265318 . 



J. I. Skullerud, P. O. Bowman, A. Kizilersu, D. B. Lein- 
weber and A. G. Wi lliams, JHEP 0304 (2003) 047 
[arXiv:hep-ph/0303176 ; A. Kizilersu, D. B. Leinweber, 
J. I. Skullerud and A. G. Wilhams, Eur. Phys. J. G 50 
(2007) 871 arXiv:hep-lat/0610078 . 
G. H. Llewellyn-Smith, Annals Phys. 53 (1969) 521. 
R. Alkofer, W. Detmold, G. S. Fischer and P. Maris , 
Phys. Rev. D 70 (2004) 014014 arXiv:hep-ph/0309077l . 
R. Alkofer, G. S. Fischer and R. Wilhams, 
arXiv:0804.3478 [hep-ph]. 

D. Zwanziger, Phys. Rev. D 65, 094039 (2002) 
[arXiv:hep-th/01 09224 . 

G. Lerche and L. von Smekal, Phys. Rev. D 65, 125006 
(2002) arXiv:hep-ph/0202194 . 

N. Nakanishi, Prog. Theor. Phys. Suppl. 43 (1969) 1. 
P. G. Tandy, Prog. Part. Nucl. Phys. 39 (1997) 117 
jarXiv:n ucl-th/9705018 . 
W. M. Yao et al. [Particle Data Group] 
(2006) 1. 

G. R. AUton, W. Armour, D. 
A. W. Thomas and R. D. Young, 
628 (2005) 125 arXiv:hep-lat/0504022; . 
D. Leinweber et al. Talk given at the conference 'Tropical 
QGD' in Port Douglas, Austraha, 27.7.-1.8.2008. 



J. Phys. G 33 

B. Leinweber, 
Phys. Lett. B 



